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Abstract 

Using the Schwinger-Keldysh closed time path integral formalism we show that the magnetic 
screening mass in QED and QCD at one loop level is non-zero as long as the single particle 

> . 

distribution function f(k) is non-isotropic, i.e. it depends on the direction of the momentum. For 
C*~) ' isotropic distribution functions such as those corresponding to thermal equilibrium the magnetic 

t>: 

screening mass at one loop level is found to be zero which is consistent with finite temperature 
field theory. The non-zero magnetic screening mass in non-isotropic non-equlibrium situations has 
Q_i| fundamental importance in that it acts as a natural cut-off to remove infrared divergences in the 



magnetic sector. Thus it allows one to avoid infrared problems which previously made it difficult 
to use a transport theory approach using perturbative QCD or QED scattering kernels to study 



the thermalization of a QED or QCD plasma. 
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I. INTRODUCTION 



An important problem to be understood in QED and QCD plasmas is how an out of 
equilibrium plasma thermalizes. This problem is quite important in the case of the QCD 
plasma produced at the RHIC and the LHC because the system lives for few fm/c after 
which it hadronizes. As the system is very short lived there is some doubt as to whether 
the system thermalizes and it is important to know the thermalization time scales of the 
various components of the plasma. As two nuclei collide almost at the speed of light the 
initial system of quarks and gluons is highly out of equilibrium. Only after many secondary 
collisions might this parton system reach local thermal equilibrium before hadronization. 
Thus there is a competition between the expansion which causes cooling and leads to the 
phase transition into hadrons, and the equilibration time scales. As all the signatures for 
quark-gluon plasma detection depend on the non-equlibrium space-time evolution of the 
partons at RHIC and LHC 0, it is essential to understand the process of the equilibration 
of partons in conditions found at these experiments. 

Due to the impracticality of solving directly using QCD the coupled Schwinger Dyson 
equations of field theory, the alternative is to solve the relativistic Boltzmann equation: 

p"dJ(x,p) = C(x,p), (1) 

in the absence of any classical background field ||, in order to study the thermalization of 
a QED or QCD plasma. In the above equation f(x,p) is the single particle distribution 
function and C(x,p) is the Boltzmann collision term describing the secondary collisions. 
The Boltzmann collision term is given by: 

C{x ^ = I wk I wk Jj0pj mpP2 pm){2 (27r)4 64{p +p *- p3 ~ ^ 

[f(x,p,)f(x,p 3 )(l ± /(x,p 2 ))(l ± f(x,p)) - f(x,p 2 )f(x,p)(l ± f{x,p A )){l ± f(x,p 3 ))]. (2) 

In the above expression M(pp 2 — > P3P4) is the matrix element for the pp 2 — > p%Pi collisions, 
(1 ± /) is due to Bose enhancement /Fermi suppression factor. Once the initial distribution 
function is known one can solve this equation to study the thermalization of the system. 
However, there is a severe infrared divergence problem that must be overcome before solv- 
ing the above transport equation. The infrared divergence comes from the matrix element 
squared \M(pp 2 — > P3P4)\ 2 when the momentum transfer: Q 2 = (p — p 3 ) 2 — > 0. This diver- 
gence occurs at small angle scattering in perturbative QED or QCD. One way to solve this 
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problem is to insert ad-hoc momentum transfer cut-off values in the integrations. However, 
the results are quite sensitive to the values of these ad-hoc cut-offs. A more satisfactory 
solution is to use a medium modified resummed propagator to evaluate the scattering ma- 
trix elements squared: \M(pp 2 —>■ P3P/t)| 2 [|J. In this way the medium dependent screening 
masses act as a natural cut-off for the infrared divergences. The calculation of the screening 
mass taking into account the contribution of all loops is an extremely difficult task and one 
therefore limits the calculation to the one or two-loop level. Although the Debye screening 
mass has been shown to be non-zero in equilibrium at one loop order, the magnetic screening 
mass is zero at one loop as shown by one loop calculations using finite temperature field 
theory methods [||, |], |7| . Hence a general supposition has been that even if the electric field 
is sceened the magnetic field is not screened at one loop level and thus one has an infrared 
problem to contend with in that approximation. 

In this paper we show that this supposition is not necessarily true. We find that the 
magnetic field is screened at the one loop level as long as the medium is non-isotropic. 
Since parton model considerations lead to non-isotropic initial conditions for the plasma, 
this leads to a natural solution to the infrared problem. All the earlier calculations done 
|5|, [7| assumed either equilibrium or isotropic distribution functions. We will show that the 
magnetic screening mass in QED and QCD at one loop level is non-zero for non-isotropic 
non-equlibrium situations. Hence there seems to be no natural infrared problem in non- 
equilibrium non-isotropic situations. This solves the problem of having infrared divergences 
in transport approaches which study the thermalization of non-equilibrium QED and QCD 
plasmas. Of course one expects that an exact study would provide non-perturbative effects 
of similar magnitude to the contribution to the magnetic mass due to non-isotropic effects. 
These effects are beyond the scope of this paper. 

The paper is organized as follows: In section II we derive the one loop magnetic and Debye 
screening mass formula in QED and QCD by using the closed-time path (CTP) formalism 
||. We summarize and conclude our main results in section III. 

II. MAGNETIC SCREENING MASS IN QED AND QCD AT ONE LOOP LEVEL 

In this section we will derive the magnetic and Debye screening mass in QED and QCD at 
the one loop level. In case of QCD we will determine only the quark loop contribution in this 
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FIG. 1: One Loop Graphs for the Gluon (Photon) Self-Energy With Quark (Electron) Loop. 

paper. The gluon loop contribution to these screening masses was discussed in a previous 
paper |I3] and we will combine these two results to discuss the total one loop contribution 
later. As we are interested in the infrared properties we will consider the infrared limit to 
define magnetic and Debye screening masses ||: i.e. we will take 



m 2 D 



ReU L ( Po = 0, \p\ -> 0) and 
m 2 q = ReU T (p Q = 0,\p\^0), 



(3) 



where U L ,U T are the longitudinal and transverse component of the self energy as defined 
below and Re means the real part of the self energy. The imaginary part of the self energy 
gives the damping rates |5], 0. 



A. Magnetic Screening Mass in QED at One Loop Level 

The diagram for the photon self energy in QED at one loop level is shown in Fig. 1. The 
retarded self energy part can be written as @: 

n£» = ^ J ^-Jr[YS u (k-p)YS 11 (k) - >fS 2X (k-p)>fS 12 (k)], (4) 

where 1, 2 correspond to +, — contour of the closed-time path. 

In this paper we will consider massless fermions since the infrared divergence problems 
occur in the massless limit. The fermion propagators S can be written as Sij(k) = ft Aij(k) 
where Aij(k) are defined below. Expressing the results in terms of retarded, advanced and 
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symmetric propagators we get from the above equation: 

n£T(p) = ^1 -0^Tr[Y^-^m^s(k- P )A R (k) + A A (k-p)A s (k) 

+ A A (k-p)A A (k) + A R (k-p)A R (k)l (5) 

where 

^ k 2 ±i sgn(ko)e 
A s (k) = -2m5(k 2 )[l - 26(k )f(k) - 29(-k )f(-k)}. (6) 

In the above equations f(k) is the non-equilibrium distribution function which must satisfy 
the constraint of yielding a finite number and energy density upon integration. As we are 
interested in studying the properties of the medium we will neglect the vacuum part of the 
self energy from the above equation. The divergence in the vacuum part of the self energy is 
absorbed in the redefinition of the bare quantities. We consider here only the terms which 
are proportional to the distribution function: we will not consider the combination A R A R 
and A A A A in eq. (^). After performing the trace in eq. (||) we get: 

IL%(p) = 2ie 2 J-0^[(k-prk v + k^k- P r-9^(k- P )-k} 

■[A s {k-p)A R (k) + A A {k - p)A s (k)}. (7) 

In the rest frame of the medium the longitudinal and transverse components of the self 
energies are given by |7|, |12| : 



Ii L {p) = L 00 (p)n 00 (p) with L o(p) = -J^; (8) 



and 



ri r (p) = 1 -T lJ (p)W(p) with r y ( P ) = -^ + p i (9) 

where we have assumed that we are interested in a static limit (po =0 and \p\ —>■ 0) so that 
we have contracted only with L 00 and in the above. If we do not consider the static limit 



we have to consider more combinations in the above equations JT2 . 

Let us first consider the transverse component IT^ (p) of the self energy: 

((* - p)V + k'(k - p)' - g"(k - p) ■ k)][A s (k - p)A R (k) + A A (k - p)A s (k)}. (10) 



Simplifying the above equation we get: 

r d 4 k 

Il£(p) = - 2.e 2 J — {% - p )k -(k.p)(p.(k- p))} 

[A s (k-p)A R (k) + A A (k-p)A s {k)). (11) 

Using the expressions for A from eq. and taking the real part we find: 

ReU T M = -2ie 2 l-^-[(k -po)k -(k-p)(p-(k-p))] 



(27T) 

o) 

-2m5(k 2 



{[ 2mS ^ P)2 \ l ~ 29(k - p )f(k -p)- 29(-k + p )f(-k + p)\ 



[l-29(k )f(k)-29(-k )f(-k)]}. (12) 



(k-p) 2 

Setting p = and again dropping the terms which does not involve the distribution function 
(the vacuum contribution) we get from the above equation: 

ReU T R (p = 0,p) = 8ne 2 J-0^[(kl-(k.p)(p.(k-p})} 

{[ 5(k l~ lk rf 2) [0(ko)f(k -p) + 9(-k )f(-k + p)] 
H - \k\ 2 

+ fT°7t lkl l^ f^ + °(-ko)f(-k))]}- (13) 
(k — \k — p] ) 

More simplification can be done by using the 5 function integration which yield: 

Aire 2 r d 3 k 



ReU T R (p = 0,p) = J 



(27r) 4 7 |fc|2_ \k-p\z 

' 'f(k-p) + f(-k + P) 



{[-\k-p\ 2 + (k-p)(p-(k-p))} 



\k-p\ 



+[\k\ 2 - (k -p)(p- (k - m l k) ) ]}• (14) 

In the limit \p\ — > we expand: f{k — p) and j^— r to the leading order and use: f{k — p) ~ 
f(k) — p ■ Vkf(k) and -J— ~ -Ml + •fi^] in the rest of our calculation. This yields: 

Attp 2 r d 3 k 1 ^ -> -> 

ReU T R {p = 0,p) ^j^-J -^ or - — {(2p- k - \p\ 2 )[f(k) + /(-*)] + 
(2tt) 4 J \k\ 2k ■ p- \p\ 2 

{-\k\ 2 - |pf + 2p ■ k + (k-p) 2 -p- k]t^[f(k) + /(-£)] - 

\k\ 2 

(1 + t^M-I^I 2 - \P\ 2 + 2p-k + (k ■ p) 2 - p- k][p- V k f(k)+p- V fe / (-£)]}. (15) 
\k\ 2 
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Taking the limit \p\ — > we get: 



m 2 (p) = ReU T R (p = 0,\p\^Q,p) 



d 3 k 



[i+(P-kn m+ J { - k) ] + ii-(p-k)Y' Vk - r p ' x " 



\k\ 



p ■ k 



p ■ k 



(16) 



After changing k — ► — k in the f(—k) part of the integration we get: 



2e k 



_fk_ 

J2n) s 



(p • 



Ifcl 



(p-fc) 2 ][ 



p ■ k 



/(*)] 



(17) 



This is the expression for the magnetic screening mass in QED at one loop level in terms of 
the non-isotropic non-equilibrium distribution function f(k). For any equilibrium distribu- 
tion function f(k) = f eq (ko) we get from the above equation m 2 g = which is consistent with 
finite temperature QED results. It can also be checked that for any isotropic distribution 
function f{\k\) we also get m 2 g = 0. Only for a non isotropic distribution function f(k) do 
we get a non-zero magnetic screening mass in QED at one loop level. 



B. Magnetic Screening Mass in QCD at One Loop Level 

First of all consider the situation of the quark loop contribution to the gluon self energy 
(see Fig. 1). The only difference between the contribution to the gluon self energy from a 
quark loop and the photon self energy with an electron loop is a factor of which comes 
from Tr[T a T b ] associated with the quark-gluon vertex. Hence the quark loop contribution 
to the magnetic mass coming from the gluon self energy is given by: 



d 3 k 

(27r)2 



\k\ 



p ■ k 



(18) 



where f q (k) is the quark distribution function. This expression is similar to the expression 



found in our earlier work [13| where we have considered the gluon loop contributions (see 



Fig. 2) in covariant gauge in the frozen ghost formalism [TCJ, |TT], |T2|. In the frozen ghost 



formalism [|10], [□], [L2|], the gluon self energy is obtained from the gluon loop and tadpole loop 
as shown in Fig. 2. In this formalism the ghost does not contribute to the medium effect 
because the initial density of states are chosen to be that of the physical gluons. It can be 
seen that apart from color factors the expression for the magnetic screening mass of quark 
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FIG. 2: One Loop Graphs for the Gluon Self-Energy with Gluon Loop and Tadpole Diagrams 

and gluon loops are the same. (Here we are only considering the contributions proportional 
to the single particle distribution function. In the vacuum sector, the quark and gluon loops 
contributions have opposite signs). 

Using Eq. fll8D and summing over all the quark flavors and then adding the gluon loop 



contributions from our previous work |13[ we obtain the final expresion for the magnetic 
screening mass in QCD: 
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d 3 k 



(2tt) 3 
+ [l-(p. kf] 



{[l + iP-kf] 



N c f g (k)+T, q f q (k) 



\k\ 



[Nj 9 (k) + T,fM 



}■ 



(19) 



p-k q 

This is the expression for the magnetic screening mass in QCD at one loop level in non- 
equilibrium. For any equilibrium distribution function f{k) = f eq (ko) we get from the above 
equation m 2 g = which is consistent with finite temperature QCD results. It can also be 
checked that for any isotropic distribution function f(\k\) we also get m 2 g = 0. Only for a 
non isotropic distribution function f(k) do we get a non-zero magnetic screening mass in 
QCD at the one loop level. 



S 



C. Debye Screening Mass in QED at One Loop Level 



Taking the (00) component of the self energy in eq. fl5|), using eqs. (||), (|3|) and carrying 
out similar algebra as above we obtain the final expression for the Debye screening mass: 

mlip) = -4e 2 / f^; Vfc /W- (20) 
J {In) 6 p . k 

This is the expression for the Debye screening mass in non-equilibrium at the one loop level 
in QED in terms of any arbitary non-isotropic distribution function f(k) of the electron. 
For any isotropic distribution function f{\k\) we get from the above equation: 

4p 2 f 

mUP) = — dk kf(k), (21) 
where k = \k\. For an equilibrium distribution function the above equation gives: 

ml = (22) 
which is the same result in QED obtained by using finite temperature QED methods || ^, [7| . 



D. Debye Screening Mass in QCD at One Loop Level 

Using the quark-gluon vertex instead of the photon-electron vertex we obtain the following 
expression for the contribution to the Debye screeing mass in QCD from a quark-loop: 

mUP) = -2, 2 / ( fi P :^ fM (23) 
J (Any p . k 

This expression for the Debye screening mass in QCD is valid in non-equilibrium situations 
for any non-isotropic distribution function f(k) for a single quark flavor. Adding the gluon 
loop expression in a covariant gauge from our earlier work Jl3| and summing over all quark 
flavors we get: 

m 2 D (p) = -2s 2 / ^^[NJM + J2 fM- (24) 
J (2vr) d p . k q 

This is our final result for the Debye screening mass in QCD at the one loop level in non- 
equilibrium for any non-isotropic distribution function f(k). For isotropic quark and gluon 
distribution functions the above equation gives: 

m%(p) = % fdkk[NJ g (k)+J2f q (k)], (25) 
9 



where k = \k\. For an equilibrium distribution of quarks and gluons at high temperature 
the above equation gives: 



ml =g 2 T 2 



No Nf 
3 6 



(26) 



which is the same result obtained in QCD by using finite temperature QCD techniques 

III. CONCLUSION 

Using the Schwinger-Keldysh closed time path integral formalism we have shown that 
the magnetic screening mass in QED and QCD at one loop level is non-zero as long as the 
distribution function is non-isotropic, i.e. it depends on the direction of the momentum. For 
isotropic distribution functions such as those describing thermal equilibrium the magnetic 
screening mass at one loop level is found to be zero which is consistent with finite temperature 
field theory results in QED and QCD. The non-zero magnetic screening mass in non-isotropic 
non-equlibrium situations has fundamental importance in that it acts as a natural cut-off 
to remove infrared divergences in the magnetic sector when we study the thermalization of 
a QED or QCD plasma using a microscopic transport theory approach. This is particularly 
important at RHIC and LHC heavy-ion collisions because it is unlikely that partons are 
thermalized quickly before hadronization. As any ad-hoc infrared cut-off put by hand change 
the plasma properties, using this natural infared cutoff we avoid ad-hoc infrared cutoff 
methods which are sensitive to the cutoff. We can use the Debye screening mass discussed 
above as the infrared cut-off in the electric sector and the magnetic mass derived in this 
paper as the infrared cut-off in the magnetic sector. Hence for non-equilibrium QED and 
QCD plasma descriptions the magnetic screening mass derived in non-equilibrium can be 
used to remove infrared divergences in the magnetic sector. 

A non-equilibrium non-perturbative calculation of the magnetic screening mass is very 
difficult. However, in ||13| , we showed that the values of the magnetic screening masses found 



by using a single particle distribtution function obtained from parton model considerations 
for conditions relevant at RHIC and LHC energies was comparable to screening masses 
obtain from lattice QCD equilibrium results at the corresponding temperatures. Hence in 
the absence of non-equilibrium non-perturbative calculations, one can use the one loop level 
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result derived in this paper to study thermalization of QED and QCD plasma using transport 
methods. 
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